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PACS 34 . 50 . -s - Scattering of atoms and molecules 
PACS 34 . 50 . Cx - Elastic; ultracold collisions 

Abstract -We use exact four-boson scattering equations in the momentum-space framework to 
study the universal properties of shallow Efimov tetramers and their dependence on the two-boson 
scattering length. We demonstrate that, in contrast to previous predictions, the shallow tetramer 
in a particular experimentally unexplored regime is not an unstable bound state but an inelastic 
virtual state. This leads to a resonant behaviour of the atom-trimer scattering length and thereby 
to a resonant enhancement of the trimer relaxation in ultracold atom-trimer mixtures. 
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Introduction. — The system of three identical bosons 
with short-range interactions in the unitary limit, i.e., 
I^when the two-particle scattering length a is infinite and 
^-Trilie dimer binding energy bd vanishes, was first studied by 
q|/. Efimov in 1970 pQ. He predicted an infinite number 
L — of zero spin and positive parity (0 + ) trimer states with 
geometric spectrum and accumulation point at the three- 
particle threshold, i.e., b n /b n —i ~ 1/515 for n — > oo, where 
f — b n is the binding energy of the nth Efimov trimer. The rea- 
^Dson for this phenomenon, called the Efimov effect, is that 
^a resonant two- particle interaction (large \a\) yields an ef- 
. fective long-range attraction in the systems with N > 3 
CO particles [T] . The resulting few-body bound states (except 
maybe for few lowest of them) are of non-clasical nature, 
|i t e.. their size is much larger than the range of the two- 
• • body interaction with the particles residing predominantly 
. m an interaction-free region. As a consequence, the prop- 
^erties of such few-body systems are universal, i.e., inde- 
$— i pendent of the details of the short-range interaction [2] . A 



characteristic feature in bosonic systems is a log-periodic 
a-dependence of few-body observables [2]. 

The experimental evidence for the Efimov physics was 
observed 35 years later in systems of ultracold atoms 
where the formation of the Efimov trimers led to a reso- 
nant enhancement or suppression of the recombination or 
relaxation processes [SHE] ; large values of a were achieved 
by variation of the external magnetic field in the vicinity of 
the Feshbach resonance. First steps in exploring the four- 
atom Efimov physics via ultracold atom experiments [SJ[7] 
have been made recently calling for accurate theoretical 



studies of the four-body systems with large \a\. 

Although there is no Efimov effect in the four-boson sys- 
tem, the three-boson Efimov effect has an impact on the 
atom-trimer scattering observablcs: scattering lengths, 
phase shifts, elastic and inelastic cross sections are related 
to the corresponding trimer binding energies in a universal 
way |S]. Furthermore, it was predicted in refs. [91ITU] that 
below each Efimov trimer two tetramers with spin/parity 
+ should exist. We label them by two integers (n, k) 
where n refers to the associated trimer and k = 1 (2) for 
a more deeply (shallowly) bound tetramer. A schematic 
representation of the four-boson energy levels as functions 
of the two-boson scattering length a is given in fig. [TJ For 
a better visualization only two families of multimers are 
shown preserving only qualitative relations between them. 
The two lowest tetramers with n = are true bound 
states; all others lie above the lowest atom-trimer thresh- 
old (ATT) and therefore are unstable bound states (UBS) 
with finite width and lifetime that could not be calculated 
in refs. [9lfT0] using standard bound state techniques. The 
widths of unstable tetramers in the unitary limit were de- 
termined using proper four-boson scattering calculations 
|S]. Trimers and the associated tetramers exist also at 
large finite \a\ as predicted in refs. [9"HTTj and confirmed in 
experiments with ultracold atoms [6j[7]. At specific large 
negative values of a where the tetramers emerge at the four 
free atom threshold they lead to a resonant enhancement 
of the four-atom recombination process [MTU] . At specific 
large positive values of a where the tetramers connect to 
the dimer-dimer threshold they lead to a resonant increase 



p-1 



A. Dcltuva 



of the dimer-dimer relaxation rate [llj . It was predicted in 
ref. [TU] that the two tetramers do not intersect the associ- 
ated ATT and lie below it in the whole region of their exis- 
tence between the four-atom and the dimer-dimer thresh- 
olds. This scenario is schematically shown in fig. [1] for 
the (n — l)th (lower) multimcr family. However, in the 
present work we prove that this is not true for the shallow 
tetramer in a particular experimentally not yet explored 
interval of positive a values where it crosses the ATT and 
becomes an inelastic virtual state (IVS); this scenario is 
schematically shown in fig. [1] for the nth (higher) multi- 
mer family. As a consequence, the atom-trimer scattering 
length shows a resonant behaviour that would lead to a 
resonant enhancement of the trimer relaxation in an ultra- 
cold mixture of atoms and excited Efimov trimers. Our 
proof relies on rigorous four-boson scattering calculations 
in the momentum-space framework. With some technical 
modifications [H] the method follows the one of refs. [T2"HT1] 
that has already been applied successfully to the descrip- 
tion of all four-nucleon elastic and transfer reactions at 
low energies. 

Atom-trimer scattering. — We use exact four- 
particle scattering equations that were introduced by Alt, 
Grassberger, and Sandhas (AGS) [H3HS]; they are equiv- 
alent to the Faddcev-Yakubovsky formulation [17] of the 
four-body problem. The AGS equations are integral equa- 
tions for the four-particle transition operators Up a ; in the 
case of identical particles they are given in refs. [5|fT2] - In 
the present work we are interested in the energy regime 
where only the atom-trimer channels are open and the 
AGS equations can formally be reduced to a single inte- 
gral equation 

U n = PuiGotGo)- 1 + U 2 (l + P 3i ) + P 3i U 1 GotG Q Un 
+U 2 G tG (l + P 34 )U 1 G tG U 11 . (1) 

The two-boson potential v enters the AGS equations via 
the two-body transition matrix t = v + vGot and the 
symmetrized 1+3 (a = 1) and 2+2 (a = 2) subsystem 
transition operators U a — P q Gq 1 + P a tGoU a - All tran- 
sition operators acquire their dependence on the available 
energy E via the free resolvent Go = (E + iQ — Hq)" 1 
with Ho being the free Hamiltonian. The permuta- 
tion operators Pjj of particles (ij) and their combinations 
Pi = P12 P23+P13 P23 and P 2 = P13 -P24 ensure the desired 
symmetry of the four-boson system. Thus, in a simplified 
interpretation the integral equation (JTJ) for the transition 
operators Up a can be viewed as a sum of the multiple 
scattering scries with all possible interactions in the four- 
particle system up to infinite order. Introducing the tran- 
sition operators t and U a at intermediate steps ensures 
that cq. ((T|) is mathematically well behaved and can be 
solved numerically Amplitudes for elastic and inelastic 
atom-trimer scattering (f\T\i) = 3(<f){\Uii\<f>\) are given 
by the on-shell matrix elements of the transition operator 
(HJ) between the Faddeev components |<^™) of the corre- 
sponding initial/final atom-trimer states. 
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Fig. 1: Schematic representation of the four-boson energy 
spectrum as a function of the two-boson scattering length. 
The tetramer (T) energies, atom-trimer (a+t), dimer-dimer 
(d+d), and dimer-atom-atom (d+a+a) thresholds are shown 
as thick solid, thin solid, dashed-dotted, and dotted curves, re- 
spectively. For a better visualization only qualitative but not 
quantitative relations between them are preserved and only 
two families of multimers, the (n — l)th and nth, are shown. 
Positive (negative) a are on the right (left) side from the ver- 
tical axis while the four free atom threshold lies on the upper 
horizontal axis. The behaviour of the (n — l,2)th tetramer 
corresponds to the predictions of ref. [10] and is non-universal 
whereas that of the (n, 2)th tetramer is universal; the dashed 
part of the curve corresponds to the IVS. 



We use the momentum-space partial-wave framework 
PEUCEB] to solve the AGS equation numerically. The 
practical solution can be simplified significantly [5] by us- 
ing an 5*-wave separable two-boson potential v = \g)X{g\. 
To prove that the universal properties of the four-boson 
system are independent of the short-range interaction de- 
tails, as in ref. [S] we performed calculations with two 
choices of the momentum-space form factor 



<p|<7) = [1 + ca (p/A) 2 ]e 



-(P/A) 2 



(2) 



namely, with c 2 — and c 2 = —9.17, that yield very dif- 
ferent off-shell behaviour of the resulting potential whose 
strength A is adjusted to reproduce the given value of a. 
Further technical details on our four-boson scattering cal- 
culations can be found in refs. [5TI12U18J . In the case of the 
true bound state (n — 0) we solve Faddeev- Yakubovsky 
equations [T7] in the version of ref. Qj|] . 

Extraction of tetramer properties. — For the clas- 
sification of the states we follow ref. [5D]. A true four- 
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boson bound state corresponds to a simple pole of the 
transition operators Up a in the physical sheet of the com- 
plex energy plane. In contrast, an UBS corresponds to 
the pole of Up a in one of the unphysical sheets adja- 
cent to the physical sheet and therefore may lead to a 
rcsonancc-likc behaviour of the scattering obscrvablcs. Let 
E n ,k = —B n ^k — *T ni fc/2 be the complex energy of the 
(n, fc)th unstable tetramer, where — B n ^ is the position 
relative to the four free particle threshold and Y n ^ is the 
width of the state. The energy dependence of Up a in the 
physical region close to E rts —B n ^ can be given by 

oo 

U Pa = J2 U { ^ k ' m \E - E n , h r (3) 

m— — 1 

where only first few terms yield nonnegligible contribu- 
tions. This allows to extract the values of B n k and T n ^ 
from the calculated atom-trimer scattering amplitudes or 
observables. Since the shallow tetramer (k = 2) is very 
close to the ATT, B n ^ and r n ^ alternatively can be ob- 
tained from the atom-trimer scattering length A n and ef- 
fective range r n [20] [21]. The elastic S- wave scattering 
amplitude as a function of the on-shell momentum k n is 
given by 

T^(k n ) = -[7r Ml (fc„cot^ - ifc„)]- 1 (4) 

where fii is the reduced atom-trimer mass. The complex 
phase shift <5„ at small momenta k n can be expanded as 
k n cot8^ = —1/A n + r n k^/2 + o(fc^). Thus, the approx- 
imate value of the complex binding momentum K n , cor- 
responding to the pole of the amplitude, can be obtained 
from the equation 

r n K 2 n /2 - iK n - 1/A n = 0. (5) 

The root K n that is closer to the threshold, i.e., smaller 
in the absolute value, relates to the tetramer position and 
width as 

- B n ,a - «r n , 2 /2 = -b n + (6) 

For the UBS ReK n < 0, lmK n > 0, and I\ 2 > [2D], 
while for the true bound state RcK n = 0, ImA„ > 0, and 
Tn.2 = 0. The method based on the effective range ex- 
pansion ([S]), although slightly less precise than the one of 
eq. (O, has an advantage of being applicable also for the 
IVS in the S wave that does not lead to a resonant be- 
haviour of the scattering observables around E ~ — -B„ j2 . 
IVS corresponds to the pole of the amplitude with 
RcA"„ < 0, ImA'„ < and r„ 2 < in the sheet of 
the complex energy plane more distant from the physi- 
cal one [20] - A particular case of the S wave IVS when 
there is no lower threshold is a virtual (unbound) state 
with RcA„ = 0, ImA„ < 0, and r„. 2 = 0. In contrast, 
the near-threshold unbound states with nonzero angular 
momentum have different properties and are called the 
Breit-Wigner resonances ^U\; a universal four-boson sys- 
tem doesn't support such states [5]. 
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Table 1: Positions and widths of shallow tetramers in the 
unitary limit. Results obtained with C2 = (c2 = —9.17) in 
eq. (|2| are given in the top (bottom) part. 

Results. — We present our results in a universal form 
such that they do not depend on the values of the particle 
mass or the cutoff parameter A. The binding energy of 
the nth excited trimer b n and the associated length scale 
L n = (2/xi6„)~ 1//2 (in our convention h = 1) will often be 
used to build dimensionlcss ratios, in particular, B n ^/b n 
and T n ^/2b n that are expected to be universal numbers. 
To achieve the universal limit we have to consider reac- 
tions with highly excited Efimov trimers that are of large 
size such that the short-range details become negligible. 
This is a serious challenge in the numerical calculations 
since the binding energies of the involved trimers differ 
by many orders of magnitude and therefore the momen- 
tum grids of correspondingly broad range must be used. 
In this work we consider reactions with up to six open 
atom-trimer channels and show that this is fully sufficient 
to obtain universal results. The predictions for the posi- 
tions and widths of shallow tetramers in the unitary limit 
are collected in table [TJ The agreement between the two 
methods based on eqs. and (|5l6p is better than 0.5% 
for n > 0. For highly excited tetramers, n > 3, the results 
approach the universal values for both choices of the form 
factor ||2J). In contrast, large deviations can be seen for 
lower tetramers where the finite range corrections become 
important. 

In the following we do not demonstrate explicitly the 
convergence of our results, however, they are checked to 
be independent of c 2 and n for n > 3 with good accuracy. 
The energy dependence of the elastic and inelastic cross 
sections for the atom scattering from the highest available 
trimer (m = n— 1) in the vicinity of the (n, 2)th tetramer, 
i.e., at E ks — -B„ j2 , is shown in fig. [2] Since the tetramer 
is UBS in the unitary limit, both cross sections exhibit a 
resonant behaviour. However, they deviate slightly from 
the exact Breit-Wigner shape due to non-negligible contri- 
butions of nonresonant (background) terms corresponding 
to m > in eq. ©. 

We consider next the shallow tetramer in the interval of 
positive two-boson scattering length a between the unitary 
limit a — > oo and the intersection of the atom-trimer and 
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Fig. 2: Elastic and inelastic cross sections for the atom-trimer 
scattering in the energy regime containing the (n, 2)th tetramer 
UBS. Unitary limit results for the highest available trimer (m = 
n — 1) are shown. 



dimcr-dimcr thresholds a = a n d where b n = 2bd- The 
evolution of the binding momentum K n with a is shown 
in fig. [3] we use the dimensionless ratio a n d /a. While 
KeK n L n remains negative and varies very slowly, ImK n 
changes the sign two times such that lmK n < for a € 

« 2 
are 



) and lmK n > elsewhere. These critical values 



a n I a n 

Thus, the shallow tetramer is an IVS for a € (a^' 2 , 



0.0729, 
0.99836. 



(7) 
(8) 



but an UBS elsewhere. The evolution of the tetramer 
properties is displayed in fig. |U Solid parts of the curves 
correspond to UBS while the dash-dotted ones to IVS. 
Starting at the unitary limit and increasing a^f/a the 
tetramer UBS moves towards the ATT while its width 
decreases. r„ i2 passes through zero at a = a n ' x where the 
tetramer becomes an IVS. By further increase of a„ /a the 
tetramer IVS moves away from the ATT but turns around 
at a re 0.65a^ d and at a = a n ' 2 where V n .2 = becomes 
an UBS again. Note that in very narrow intervals close 
to a = a n ^ the tetramer (in both UBS and IVS cases) 
is slightly above the ATT, i.e., B„,2 < b n . This is evi- 
dent from eq. © with Imi"^ = and can be seen in the 
topinset of fig. |U 

We do not show the evolution of the tetramer properties 
for aff Ja > 1 where our preliminary results are in quali- 
tative agreement with ref. |llj : the (n, 2)th tetramer is an 
UBS with B n 2 > 2bd and finite width until it intersects 
the dimcr-dimcr threshold at a 
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The tetramer IVS resides on a sheet of the complex 
energy plane that is not adjacent to the physical one. It 
therefore affects the atom-trimer scattering observables in 
a completely different way as compared to UBS: the cross 
sections around E « —B n 2 do not exhibit the resonance- 
like behaviour seen in fig. [H but may have a cusp at the 
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Fig. 3: Real and imaginary parts of the complex binding mo- 
mentum K n of the shallow tetramer as functions of the two- 
boson scattering length a. 
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Fig. 4: Position of the shallow tetramer relative to the ATT 
(top) and its width (bottom) as functions of the two-boson 
scattering length a. The thin dotted line represents the zero 
value. 



ATT (E = -b n ) if the IVS pole is very close to it. This is 
illustrated in fig. [S] where the atom-trimer inelastic (n — 
1 — > n — 2) cross section is shown for two cases. At /a 
0.0820 the tetramer IVS with (b n - B na - iT n>2 /2)/b n w 
(-1.93 + i2.44) x 10" 5 is very close to the ATT while 
it is more distant at a n d /a ~ 0.1682 with {b n — B n 2 — 
ir„ >2 /2)/6„ re (-2.04 + i0.22) x 10~ 3 . As a consequence, 
the cusp at the ATT is clearly pronounced in the first 
case but hardly visible in the second. We emphasize that 
such behaviour is characteristic to S waves only as the 
unbound states with nonzero angular momentum close to 
the threshold are the Breit-Wigner resonances. 

The resonant peak for the UBS and the cusp for the 
IVS can be understood given the pole factor in the ampli- 
tude (U) that is common for all elastic and inelastic chan- 
nels. For simplicity we consider the case \k n \ 1 \K n \ « 
l/|r„|, i.e., when the available scattering energy and the 
UBS/IVS poles are very close to the nth ATT. Then the 
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Fig. 5: Cross section for the inelastic atom-trimer (m = (n — 
l)th state) scattering leading to the (m — l)th trimer state in 
the vicinity of the nth ATT. Two cases differing in the complex 
energy of the (n, 2)th tetramer IVS are shown. 
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Fig. 6: Atom-trimer scattering length A n as function of the 
two-boson scattering length a. 



energy-dependence of the cross sections (up to channel- 
dependent factors) is approximately given by 

cr - \K n -k n \~ 2 

= [(RcK n ) 2 + (lmK n - y/2m\b n + El) 2 }- 1 

x 9(-6„ - E) 

+ \{lmK n f + (ReK n - + El) 2 }- 1 

x Q(b n + E). (9) 

Here Q(x) is the Heaviside step function with values 1 and 
for x > and x < 0, respectively. Given ReK n < it 
is easy to see that a resonant peak only takes place for 
ImK n > at E < —b n , i.e., for the UBS below the nth 
ATT. In contrast, for lmK n < the cross sections reach 
maximum value at E = — b n that corresponds to the IVS 
cusp at the ATT. 

The two intersections of the ATT with the tetramer 
at a = a^ J lead to a resonant behaviour of the atom- 
trimer scattering length A n as demonstrated in fig. El The 
positive (negative) peaks of RcA n correspond to tetramer 
UBS (IVS) with B nt 2 = b„ while ReA n = and negative 
peaks of lm.A n correspond to T„ 2 = 0. No experiments 
have been performed so far in this regime but the UBS- 
IVS conversion could be observed by creating an ultracold 
mixture of atoms and excited Efimov trimers in a trap 
and measuring the trimer relaxation as function of a. Ac- 
cording to ref. [5] , the zero-temperature limit of the trimer 
relaxation rate constant is /3° = — (47r//ii) JmA n . Thus, 
a resonant enhancement of the trimer relaxation process 
takes place at a — that should be visible at sufficiently 
low temperature, possibly partially smeared out due to the 
thermal averaging. On the other hand, in the real exper- 
iments deviations of from the universal values ([7][5]) 
can be expected due to finite-range effects and/or presence 
of deeply bound dimers. 



Comparison with previous works. Finally, we 
would like to discuss the differences between our results 
and those of ref. [TD] . At the unitary limit ref. [TU] pre- 
dicts {Bn^/bn — 1) ~ 0.01 to be more than 4 times larger 
than our universal result 0.00228 presented in tableQ] The 
results of ref. [10] were obtained using the hyperspheri- 
cal framework in coordinate space and neglecting the fi- 
nite width of the tetramcr. In addition, for the (n, 2)th 
tetramer they were limited to n < 1 (n = with strong 
repulsive three-body force) and by far have not reached 
the level of convergence comparable to ours. For ex- 
ample, our five best-converged results for (B n ^/b n — 1) 
are 0.00228 within 0.5% accuracy while the three best- 
converged results of ref. [TU] are 0.006, 0.03, and 0.001. 
On the other hand, our n < 1 results given in table [JJ 
that receive nonnegligible finite-range corrections show de- 
viations from the universal limit comparable to those of 
ref. [TU]. In the presence of finite-range effects we found 
that increasing (B n ^/b n — 1) at a —> oo shifts a^f/a^' 1 
to larger values thereby shrinking the IVS region. Thus, 
when the tetramcr at a — ¥ oo is sufficiently far from the 
ATT, it may even remain UBS (true bound state if n = 0) 
in the whole regime a^ja < 1. Such a situation takes 
place also for the non-universal (0, 2)th tetramer of the 
present work. However, by choosing ci > 1 in eq. (J2J 
we obtain the (0, 2)th tetramer that mimics the universal 
behaviour, i.e., makes transition into a virtual state. This 
proves again that low n states are sensitive to the short- 
range details and may explain why the tetramer IVS could 
not be predicted in ref. [TU] . 

In the case of the deeper (k = 1) tetramer there is a rea- 
sonanble agreement between our value B n ,i/b n = 4.611(1) 
obtained in ref. [5] and the one of ref. [TU], 4.58, the latter 
having about 1% uncertainty. 
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Summary. — Wc studied the universal properties of 
the shallow Efimov tctramcr. Since this four-boson state 
lies in the continuum, we solved exact AGS four-particle 
equations for the atom-trimer scattering. The momentum- 
space partial-wave framework was employed. The univer- 
sal limit of the results is achieved accurately in reactions 
with highly excited trimers (up to 5th excited state). Po- 
sitions and widths of shallow tetramers were calculated as 
functions of the two-boson scattering length. We demon- 
strated that, in contrast to previous predictions by other 
authors 01101 obtained with some limitations, the uni- 
versal shallow tetramer intersects the ATT twice and in 
a regime between those intersections a e (aJJ; 2 , a^' 1 ) it is 
an inelastic virtual state. The tetramer is an unstable 
bound state outside the above interval of a, in qualitative 
agreement with the results of refs. [TUHH]. We studied 
how the atom-trimer scattering obscrvables are affected 
by these changes of the tetramer properties. In particu- 
lar, wc demonstrated a resonant behaviour of the atom- 
trimer scattering length that could be observed as a reso- 
nant enhancement of the trimer relaxation in the ultracold 
mixture of atoms and excited trimers. 

* * * 

The author thanks R. Lazauskas for discussions and 
A. C. Fonseca for comments on the manuscript. 
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